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Supersymmetric models with spontaneously broken approximate R-symmetry contain a light spin 
particle, the R-axion. The properties of the particle can be a powerful probe of the structure of 
the new physics. In this paper, we discuss the possibilities of the R-axion detection at the LHC 
experiments. It is challenge to observe this light particle in the LHC environment. However, for 
typical values in which the mass of the R-axion is a few hundred MeV, we show that those particles 
can be detected by searching for displaced vertices from R-axion decay. 



I. INTRODUCTION 



i— h . Discovering supersymmetry is one of the main goal of the LHC experiments. Once the supersymmetric standard 
I ^ ' model (SSM) is confirmed experimentally, the next question is how the supersymmetry is broken and how the effects 
of symmetry breaking are mediated to the SSM sector. In most cases, such investigations on "beyond the SSM 
physics" rely on arguments based on extrapolations of the observed supersymmetry mass parameters to higher energies, 
i— - n , However, there is one class of models of supersymmetry breaking in which we can get a glimpse of the structure of 

• the hidden sector in a direct way with the help of the R-symmetry. 
Oh| The R-symmetry is a symmetry which is rather generic in the models of spontaneous supersymmetry breaking 
through _F-term[l]. At the same time, however, it must be broken in some way in order for the gauginos in the 
SSM sector to have non-vanishing masses, which we are assuming to be the Majorana fermions. One possibility of 
, H the gaugino mass generation is to consider models where the gaugino masses are generated as a result of the explicit 
breaking of the R-symmetries, which is realized in the supergravity mediation scenario [2, 3], the anomaly mediation 
, scenarios [4], and in a class of the gauge mediation [5-8] with the explicit R-symmetry breaking [9-14]. Unfortunately, 
■ in those models, the R-symmetry leaves little trace for the collider experiments, since the mass of the R-axion is 
CO ' typically heavy and beyond the reach of the LHC experiments. 

In this paper, instead, we focus on a class of models with gauge mediation where the R-symmetry in the hid- 
den/messenger sectors is exact in the limit of the infinite reduced Planck scale, i.e. Mpl — > oo. Then, the gaugino 
masses are generated as a result of the spontaneous breaking of the R-symmetry. We also assume that the R-symmetry 
is respected by the SSM sector as well as the sector of the origin of the higgsino mass \i and the Higgs mass mixing 
B[i except for the anomaly to the SSM gauge interactions. We call this scenario, the minimal R-symmetry breaking 
scenario. 

The immediate consequence of the spontaneously broken approximate R-symmetry is the existence of a pseudo 
Nambu-Goldstone boson, so called the R-axion [1]. As is well know, the R-axion can have a large mass [15] compared 
k^l ' with the Peccei-Quinn axion [16, 17], while getting small couplings. Hence, it is free from cosmological and astrophys- 
ical constraints. The notable feature of the R-axion in the minimal R-symmetry breaking scenario is that it couples 
to the gluinos and this enhances the R-axion-gluon coupling at the low energy effective theory. As we will show, a 
sizable number of R-axions can be produced via the gluon fusion process with the dominant contribution coming from 
gluino and messenger loops. Since the R-axion has long lifetime, it is possible to detect the R-axion by searching for 
the displaced vertex of the R-axion decay. 

It is possible that the R-axion can produce the following striking signature: A promptly produced jet of px > 
100 GeV, recoiling against a low mass muon pair produced in the beam pipe vacuum region. In the paper, we will 
compute the rate of background for this signature. 

The organization of the paper is as follows. In the next section, we briefly review the generic properties of R- 
symmetry breaking in models with gauge mediation. In section III, we summarize interactions of the R-axion with 
the SSM particles. Now the paper shifts from a theoretical to an experimental viewpoint. In section IV, we discuss 
the production of the R-axion and the typical signatures. In section V, we give a rough estimation of the background 
processes to the R-axion signal. The final section gives our conclusions. 
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II. R-SYMMETRY IN GAUGE MEDIATION SCENARIO 

As examples of models of the R-symmetry breaking, we consider two simple examples of the models with gauge 
mediation. The first example is the so-called 'minimal gauge mediation' model in which the supersymmetry breaking 
effects are mediated by a pair messenger fields ip and f/> which are charged under the SSM gauge group. To handle the 
supersymmetry breaking, we introduce a chiral superfield S whose F-term has a non-vanishing vacuum expectation 
value (VEV) called F in the following discussion. In the minimal model, the messenger fields and the supersymmetry 
breaking field couple through the superpotential, 

W = kSipip, (1) 

where k is a coupling constant. The above superpotential possesses an R-symmetry with a charge assignment S(2), 
tp(0), and ip(0), which forbids the gaugino mass if it is not broken. This charge assignment is consistent with F ^ 0, 
since it implies that the F-term has R-charge 0; other possible ?7(l)'s, consistent with nonzero gaugino masses, are 
broken by F. 

In this example, the scalar VEV of the chiral superfield S provides a supersymmetric mass term to the messenger 
particles, which corresponds to the order parameter of the R-symmetry. Then, in order that the messenger fields are 
not tachyonic, the scalar and the F-term VEVs of the chiral superfield S, 

(S)=M + F6 2 , (2) 

are required to satisfy the condition, 

\kM\ > \y/F\. (3) 

Thus, in the minimal gauge mediation, the R-symmetry must be broken, with a VEV even larger than the supersym- 
metry breaking scale as long as k is perturbative. 

Another example is a model proposed in Ref. [18]. Let us add another pair of the messenger particles tp' and xp' so 
that we can give the supersymmetric and R-symmetric mass terms to the original messengers; 

W = kSipip + rmpft + mip'-ip, (4) 

where m denotes a mass parameter. Here, again, the R-symmetry which is relevant for the gaugino mass suppression 
is defined by S(2), -0(0), V>(0), and 4>'{2). In this model, the messenger fields are not tachyonic even if M = 

contrary to the former example. However, even in this case, there is a lower bound on the size of M coming from the 
experimental lower bound on the gaugino masses. That is, in this model, the gaugino mass is roughly given by, 



a kM 
47r m 



kF 



m 

,2 



kF* 

(5) 



Then, remembering that \kF/m 2 \ < 1 so that the messengers are not tachyonic, we obtain a lower bound on M, 



4-7T 

M > — m gaugino , (6) 



where the lower bound is saturated for m ~ a/F- 

The above discussion applies to general gauge mediation models with R-symmetric hidden/messenger sectors. In a 
general gauge mediation model, the gaugino masses can be parameterized by, 

a M F* ( ( F \\ 

™~ * Nm ^M^MZ [ Cm + ° \Ml ) ) ' (?) 

where the coefficient C m is a model dependent parameter at most O(l), M m (> M) is the scale of the gauge mediation, 
and N m is the number of the messenger fields. 1 Here, we have omitted 0(1) coupling constants of the mediation 
mechanism. From this expression, we obtain a lower bound on the R-symmetry breaking scale, 

tfgtoW'VJSgL). (8) 



1 For recent developments on the analysis of the models with the general gauge mediation, see Refs. [19—23]. 
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which results from the non-tachyonic messenger condition, M m > \[F ? The lower bound is saturated for the lowest 
scale gauge mediation, i.e. M m ~ \/~F. Therefore, we see that the R-symmetry breaking scale is about 10 4 GeV or 
higher in generic gauge mediation models. 

Notice that the masses of the SSM scalars are roughly given by, 

I a F 

^scalar ~ V N m — — — , (9) 
47T M m 

and that there is no suppression from the R-symmetry breaking. Thus, when the R-symmetry breaking scale is 
suppressed compared with M m and \/F, the gaugino masses are suppressed compared with the scalar masses. In 
other words, the SSM scalars are much heavier than the electroweak scale even if the gauginos have the masses of the 
order of the electroweak scale. However, such heavy scalars bring up the hierarchy problem again; we do not pursue 
the possibility of such a highly hierarchical spectrum in this study. Instead, we consider models where the gaugino 
and scalar masses are not so different, i.e. M ~ M m ~ y/F, which is naturally realized in low scale direct gauge 
mediation models (see for example [24-27]). 

In the following, we confine ourselves to the scenario with the low R-symmetry breaking scale as well as the low 
scale supersymmetry breaking and the messenger scale, i.e the lowest gauge mediation scenario; 

fn = 10 4 - 5 GeV, 
Vf = 10 4 - 5 GeV, 

M m = 10 4 - 5 GeV, (10) 

where Jr denotes the R-symmetry breaking scale. As we will show in the subsequent sections, a sizable number of 
R-axions can be produced in the LHC experiments for such low scale gauge mediation models. 

We mention here that the low scale gauge mediation scenario is also motivated from cosmology. The low scale 
gauge mediation scenario given above admits a very light gravitino with the mass lighter than O(10) eV. Such a light 
gravitino is free from all the cosmological problems [28] . The model has the phenomenology typical in gauge mediation 
that the SSM particle is unstable and decays to the gravitino. 



R-axion Mass 



One immediate consequence of the spontaneous breaking of the R-symmetry breaking is the R-axion. In our 
scenario, we have assumed that, if supergravity in the action is ignored (Mpl — > oo), the R-symmctry is only broken 
by the anomalies of the SSM gauge interactions. In this limit, the R-axion obtains a mass only from the QCD 
anomaly. However, supergravity adds several R-symmetry breaking effects which contribute to the mass of the R- 
axion [1] . In particular, the R-symmetry is broken by the constant term in the superpotential which is always required 
to set the cosmological constant to zero after the supersymmetry breaking. 3 For example, in the case of the simplest 
supersymmetry breaking sector model with the full superpotential, 

W = A 2 susy S + w, (11) 

the R-symmctry is explicitly broken by the constant term w. Here A susy denotes a dimensionful constant such that 
the F-term VEV of S is given by (F) = A 2 usy , and the flat universe condition requires A 4 usy = 3w 2 /Mp L . In the 
supergravity, this superpotential leads to the R-symmetry breaking term in the scalar potential of S, 

V i = - 2 Ml k —y S + C - C ' = - 2 ™3/ 2 A s 2 usy S + c.c. (12) 

Here we have used the definition of the gravitino mass m 3 / 2 = w/M p . S obtains a VEV, and we parametrize this by, 

S = fR/y/Se*"*, f R then gives the scale of spontaneous R-symmetry breaking. As a result of the breaking term 
V^, the R-axion obtains a non- vanishing mass 

2 2V2m 3/2 A s 2 usy 2V6ml M PL 

™> a = = • (13) 

JR JR 



2 The number of the messenger fields is consistent with the perturbative grand unification theory (GUT) as long as it satisfies a condition, 
N m < f50/ln(M GUT /M 

mess ) • 

3 We may further postulate that the constant term also emerges as a result of the spontaneous symmetry breaking of another R-symmctry. 
In this case, we have an additional axion. The additional axion can provide the solution to the strong CP-problem. 
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Thus, if we consider the low breaking scale and the light gravitino mass, the R-axion mass is given by, 

—« m « v (^)(to4v)" 2 ' <»» 

In general models, we could consider other explicit R-symmetry breaking terms that are also vanishing in the limit 
(see Ref. [15] for more discussion). Those breaking terms also contribute to the R-axion mass at the similar size of the 
above contribution especially for the lowest gauge mediation scenario. We should also notice that, in general models, 
there is a distribution between the total supersymmetry breaking scale, A susy , which determines the gravitino mass, 
and the scale used in the gauge meditation, \/F. Therefore, the R-axion mass depends on the details of the models 
of the hidden/messenger sectors and could be larger than the value given in Eq. (14). 4 

In the following analysis, we take the R-axion as an independent parameter. We mainly focus on the range of 
R-axion masses, 

m a = 0(100) MeV, (15) 

which is, for example, the most likely mass range for f R <~ y/F > 10 4 GeV and a light gravitino mass, to 3 / 2 = O(10) eV 
in the above simple supersymmetry breaking scenario. As we will see in the next section, R-axions in this range have 
long lifetimes and leave displaced vertices inside the detecters. 



III. INTERACTIONS OF R-AXION 



Let us now consider the interactions between the R-axion and the SSM sector. Here, we consider the minimal Higgs 
sector model with two Higgs doublets (i.e. the MSSM), although we can easily extend our analysis to models with 
a non-minimal Higgs sector. As we will see, the R-axion mixes with the pseudoscalar components in the two Higgs 
doublets at the energy scale below the electroweak symmetry breaking. 5 The mixing angles determine the couplings 
between the R-axion and the SSM matter fields. 

For our purpose, the effective Lagrangian approach valid below the messenger scale is the most useful. After 
integrating out the hidden/messenger sectors, the low energy effective theory is given by the SSM particles and the 
R-axion a. In the effective theory, the R-symmetry is realized by the shift symmetry of the R-axion, 

a/fa - a'/f R = a/f R + 2c, (16) 

with the angle c, as well as by the linear transformation of the SSM fields with the charges given in Table I. Notice 
that the R-charges of the SSM fields satisfy, 

X u + X d = 2, X u + Xq + Xjj = 2, X d + Xq + Xq = 2, (17) 

so that the R-symmetry is consistent with the /i-term and the Yukawa interactions. 6 

In the effective theory just below the messenger scale, the R-axion only emerges in the gaugino mass terms, and 
the anomaly coupling terms. The effective Lagrangian of the R-axion is given by, 

£~a = \{d~af - m\ /2 e-^f«W + C H ^^-F l F* + c.c, (18) 

where gi denotes the gauge coupling constant, m\^ 2 the gaugino mass, F^ v = e^ vpa F / " 7 /2, and i runs over the SSM 

gauge groups. The coefficient of the anomaly coupling terms Cr depends on the model of the hidden/messenger 
sectors. From the second term, we find that the R-axion couples to the gauginos by a Yukawa-type interaction, 

A„t = -i-^-aXX + c.c. (19) 

JR 



4 In models in [18, 29], the R-axion can be much lighter than the values discussed above, thanks to the scale invariance of the superpotcntial 
of the supersymmetry breaking field. See also Refs. [30] for related works. 

5 In this study, we assume that there is no supersymmctric CP-problem, which is one of the motivation of the models with the gauge 
mediation. 

6 In this study, we assume that the /i-term carries a vanishing R-charge. In models where the /i-term is also generated via the interactions 
between the hidden/messenger sectors and the Higgs sector, however, the resultant /t-term can have non- vanishing R-charge. The 
following analysis can be generalized straightforwardly to such models. 
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TABLE I: The R-charges of the supermultiplets in the SSM and gauginos. Here, we also showed the charges of them under the 
Standard Model gauge group; SU(2) xU(l)y. The charge assignment is identical for all the three generations. 



At a lower energy scale, the R-axion dependence in the gaugino masses spreads to the other soft supersymmetry 
breaking masses, such as the B/i and the A-terms, as a result of renormalization group evolution. The R-axion 
dependence in the B/i-term plays a particularly important role in determining the mixing between the R-axion and 
the SSM Higgs bosons. The relevant terms for the R-axion-Higgs mixing are given by, 

V = (H 2 +m 2 H J\Hl\ + {^\ 2 +m 2 Hd )\Hl\ - {eT^B^H^ + c.c.) + \{g 2 + g' 2 ){\H° u \ 2 - \H°\ 2 ) 2 , (20) 

where H® d are the QED neutral components in the Higgs doublets. Notice that in these models, the B\i term is 
generated by the renormalization group evolution coming from gaugino loop diagrams [31] and so is proportional to 
e -ia/f R j rp Q ana jy ze the mixing mass matrix, it is convenient to decompose the neutral Higgs fields into the radial 
and the axial components, 

Hi = -±=K + p u )e^, H° d = ^=(v d + p d y^. (21) 

Here, v u ,d are vacuum expectation values of the Higgs fields which are related to the Z° boson mass by 

v 2 ee vl + vj = 4m%/(g 2 + g' 2 ) ~ (246 GeV) 2 . (22) 
The ratio between v u and Vd is traditionally expressed by, 

tan 13 ee v u /v d . (23) 
Then, with the aid of the minimum conditions of the and H®; 

\lA 2 +m 2 Hu = B/j,cot(3+ (m|/2)cos2/3, 

\fi\ 2 +m 2 H ] = B^tan/3- (m|/2) cos 2/3, (24) 

we obtain the mixing mass matrix of the axial components and the R-axion, 

i f L\ ( cot/3 1 -rcos/3 \ 

V mix = -x*M 2 x, x= \ £,d , M 2 = Bn\ 1 tan/3 -rsin/3 , (25) 

\a J \ —rcos/3 —rsin/3 r 2 cos/3sin/3 / 

where r is defined by r = w//fl- 8 

As we expected, the above mass matrix possesses two massless states, i.e. rank(Al 2 ) = 1. One of them corresponds 
to the would-be Goldstone boson absorbed by the U(1)y gauge boson, and the other one is the R-axion. The mass 
eigenstates of the mass matrix are given by, 

sin/3 -cos/3 \ / & \ 

kcos/3 Ksin/3 — /-cr sin /3 cos /3 II £ d I (26) 

nr cos 2 /3 sin /3 nr sin 2 /3 cos /3 k J y a J 

(m 2 a ,m\ ,ml) =(o,J§^,o), (27) 







? 


- ( 







The Bfi term also gets small contribution from the two-loop diagrams in which the messengers circulate as well as the gauginos and 
Higgsinos [31], which is also proportional to e~ la ^ R . 

In this study, we fix B/j, > which is realized by using the definitions of R-symmctry and the Pcccci-Quinn symmetry. 
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where k is defined by k = (1 + r 2 sin 2 2/3) x l 2 . As a result, the axial parts of the Higgs fields mix with the low energy 
R-axion; 

~ nr cos 2 (3 sin (3 x a, £<j ~ nr sin 2 /3 cos (3 x a. (28) 

The interactions of the R-axion with the SM fermions are the same found in the DFS axion model [32] . From the 
mixing in Eq. (28), the coupling constants between the R-axion and the SM fermions are given by, 

K = iyu/V2r cos 2 (3 sin (3 = im u /f R cos 2 (3, X d = iy d /V2r sin 2 (3 cos (3 = im d /f R sin 2 (3, (29) 

for up- and down- type quarks (and for their higher generation counterparts), respectively. The leptons also couple 
to the R-axion in a similar way of the down- type quarks. Unlike the case of the CP-odd Higgs scalar Ao, there is no 
tan/3 enhancement in coupling constants of the down-type quarks. 

Notice that, for a given R-charge assignment to the /x-term, the R-charges of the Higgs bosons are determined 
unambiguously so that the would-be-goldstone boson Go is invariant under the R-symmetry. That is, under the 
U(l) R symmetry, the axial fields £ Mj< 2 are sifted by, 

£u/v u ^ in jv u = i u lv u + X u c, £ d / v d — > (, d '/v d = £d/v d + X d c. (30) 

Then, by requiring that the would-be goldstone boson Go is invariant under the R-symmetry, we obtain a condition, 

X„sin 2 /3-X d cos 2 /3 = 0. 

By remembering the conditions in Eq. (17), we find that the R-charges are given by, 

X u = 2 cos 2 (3, X d = 2 sin 2 (3. (31) 

We can also check that A is also invariant under the R-symmetry with the above charge assignment, and only the 
R-axion a in the low energy theory is shifted by, 

a -> a' = a + 2«T 1 /flC. (32) 

This means that the decay constant of the R-axion below the electroweak scale is changed from f R to 

!a = K~ l f R > f R , (33) 

although k — 1 and /a — f R for f R 3> v. 

Finally, let us comment on the coefficient of the anomaly coupling in the low energy effective theory. As we 
mentioned above, the coefficient Ch depends on the model of the messenger sectors. For example, in the case of 
minimal gauge mediation with N m messenger pairs, the coefficient is given by Ch = —N m , while Ch — in the case 
of the second model discussed in the previous section. 9 At lower energy, the coefficient gets additional contributions 
from heavy fermions such as gauginos and top quarks corresponding to the R-charges of their species. Thus, at the 
scale below the gaugino masses and the top-quark mass, the coefficient of the anomaly term is changed to 

£ cff = C L ^^FF = (C H + G 2 (G) - cos 2 p)J?J?-FF, (34) 

which reproduces the anomaly of the R-symmetry by the shift of a. Here Cz(G) denotes the quadratic Casimir 
invariant which corresponds to the gauginos, and we have used the charges of the Higgs bosons fixed in Eq. (31) to 
determine the contribution from the top quarks in the third term. 

At a lower energy scale, the lighter fermions also contribute to the anomaly coupling once they are integrated out. 
For example, for the R-axion with a hundred-MeV mass, the value of the anomaly coefficient to photons below all 
hadron thresholds, is given by, 

G L , 7 = C H ,i + C h ,2 + 2-9 {^X u + 1 X^j -X d -12 tr (Q A Q 2 em ) , (35) 



9 The models proposed in Refs. [26, 33, 34] also give a vanishing coefficient, Ch- 
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where the first two terms denote the messenger contributions to the U(1)y and SU(2)l gauge bosons, the third term 
comes from the SU(2)l gaugino contribution, and the fourth from the tau lepton contribution. The final term comes 
from the light quark rotation with a charge assignment Qa which eliminates the anomalous coupling to the gluons [35] , 
and the explicit charge assignment of this manipulation is given by, 

Qa = ^|^M- 1 /tr(M- 1 ) ~ 0.16 C L , gluon , (36) 

where M q denotes the quark mass matrix and we have used the mass ratio, m u /md ~ 0.56 in the last expression. The 
coefficient Cl, giuon denotes the anomalous coupling to the gluons at the low energy scale, 

3 

Cl, giuon = Ch + 3 — - (X u + X d ) . (37) 

Therefore, we obtain those coefficients, 

Cl, 7 — Ch.i + Ch,2 — 2 — 1.9 Ch,3, 
Cl, giuon — Ch,3, (38) 

for a large value of tan/3. 



Decay properties of R-axion 

In this subsection, we analyze the decay properties of the R-axion with a hundred-MeV mass based on the effective 
Lagrangian obtained above. In this mass region, the available decay modes of the R-axion are those into a muon pair, 
an electron pair, a photon pair, and a gravitino pair. As we will show below, the dominant decay modes are the one 
into the muons for m a > 2m ^, while the one into an electron pair dominates for m a < 2m ^. 

The R-axion decays into a pair of the SM fermions via the Yukawa interaction given in Eq. (29) and the decay 
width is given by, 




Thus, the decay rate and the decay length at the rest frame of the electron pair mode are given by 
T ee ^ 3.1 x 10~ 17 GeV x sin 4 /? (- 



. OS 1/2 

4m; x ' 



,300McV 



^^■^lO^x^Wl-^-l 11-^1 . 



sin 4 /? 

while the ones into a pair of the muons are given by 



r w ~ 1.3 x 10~ 12 GeV x sin 4 /3 ( 




= « =< ""-cm x ^jjfg^l | 1 - | ■ (41) 

The R-axion also decays into a pair of photons via the anomaly coupling as appeared in Eq. (34). The width of this 
mode is given by, 



C 



2 



16tt 



f — Y ( T^L 

UttJ \f R 



r 77 ~ -^p- ( — ) [ ^ ) m. 



— "»-M^)'(T)- < 42 » 

where, Cl i7 is given in Eq. (38). Compared with the decay widths into a pair of the electrons and muons given above, 
this process is subdominant for m a = 0(100) MeV. 
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Since we are considering the light gravitino scenario, the R-axion also decays into a gravitino pair, ipip, with the 
decay rate [15], 10 

r 1 fm 3/2 \ 2 

r# - [it ) ma ' 

As a result, the decay mode into a pair of the gravitinos is also subdominant compared with the modes into electrons 
and muons, as long as we are considering the light gravitino scenario, m.3/2 = O(10) eV. 

Before closing this section, we comment on the decay properties of the R-axion with much heavier mass. When 
the R-axion mass gets heavier than 3m T , the decay mode into three pions is open. Since the direct coupling between 
the R-axion and the u, <i-quarks are suppressed (see Eq. (29)), the dominant contribution to the decay width of the 
three pion mode comes from the anomaly coupling; aFF. However, we found that the decay width into three pions 
is always suppressed compared with the one into a muon pair for m a < 1 GeV by using an effective Lagrangian of the 
R-axion— 3n interaction given in Ref. [36]. 

For much heavier R-axion, the decay modes into heavier fermions, such as the r lepton and the bottom quark 
become available. In such a region, the decay length is much shorter than the one discussed above. Thus, the R-axion 
detection in the heavier mass region is more challenging at the LHC than the case of the lighter R-axion with a mass 
in the hundred MeV range. 



IV. R-AXION PRODUCTION AT THE LHC 



In this section, we consider R-axion production with a large transverse momentum at the LHC via the gluon fusion 
processes gg — > ga, qg — > qa, and qq — > ga. These processes are similar to the CP-odd Higgs boson production 
processes in the MSSM [37, 38]. In those processes, the R-axion couples to the gluons via the anomaly coupling given 
in Eq. (18) as well as the triangle and box diagrams in which the gluino and the quarks circulate. When the transverse 
momentum pt of the R-axion is smaller than the gluino mass the anomaly couplings in the effective Lagrangian in 
Eq. (34) give a good approximation to the gluino-loop diagrams. Notice that the top-loop contribution is proportional 
to cos 2 /?, and hence, the top-quark contribution is suppressed for tan (3 > 2 to the gluon fusion processes. As for the 
contributions from the other lighter quarks, they are also suppressed compared with the gluino contribution when pt 
is much larger than the masses of those lighter quarks [39] . 

From the above discussion, we consider the R-axion production via the effective anomaly coupling, 

£ - = Cl Sj-u f ^ v = + N ^Si F ^ (44) 

where N c — 3 is a color factor. This should give a good approximation for pt < m g i u i no . The matrix elements of the 
processes, gg — > ga, qg — ► qa, and qq — > ga, from the above effective anomaly coupling are given in Refs. [37, 38]. 

In Fig. 1, we present the transverse momentum distribution of the R-axion event with the mono-jet, dan/dpr at 
the LHC and the Tevatron, for fx = 10 4 GeV, and m a — 300 MeV. In our analysis, we have used the CTEQ5 parton 
distribution functions [40] with the QCD scale Ag = 165.2 MeV. The factorization and renormalization scales have 
been set to px- The model dependent parameter Cl has been set to C'l = 2 which corresponds to the minimal gauge 
mediation model with N m — 1 and N m — 5 (up to the sign of Cl). The transverse momentum distribution is almost 
independent of the R-axion mass as long as m a < 1 GeV, while it scales with f^ 2 and C\. The figure shows that a 
sizable number of the R-axion events at the LHC even for pr > 100 GeV, while the cross section is rather suppressed 
at the Tevatron. 

Now the question is how to detect the R-axion events. As we have seen in the previous section, the R-axion has 
a long lifetime. A useful variable is transverse decay length Lt — -fcrsinO, where 7 and 9 are the boost factor 
and the polar angle of the R-axion, respectively. To see how the R-axion behaves after its production, we have 
plotted the transverse decay length distribution of the R-axion decaying into a pair of muons for m a = 300 MeV in 
Fig. 1. The figure shows that the produced R-axion mainly decays at around O(l) cm, O(10) cm, and O(100) cm, for 



For example, the helicity unsuppressed operators similar to the gaugino coupling in Eq. (18) is proportional to e la ^Ripip at the leading 
order, which is forbidden in the masslcss gravitino and the exact R-symmetry limit. 
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FIG. 1: Left) Transverse momentum distribution of the R-axion with the mono-jet event at the LHC (solid) and the Tevatron 
(dashed). Here, we have taken }r = 10 4 GeV , m a = 300 MeV, and Cl = 2. We have integrated the rapidity y for \y\ < 3. 
Right) Transverse decay length distribution of the R-axion decaying into a pair of the muons at the LHC. The solid lines denote 
the distributions for the minimum transverse momentum, pr.min = 100 GeV for fu = 10 4,4 5,5 GeV from the bottom up. The 
dashed lines correspond to pT,min = 200 GeV, and the dotted line to pr,min = 300 GeV for fa = 10 4 ' 4,5,5 GeV from the bottom 
up. 



Jr = 10 4 GeV, Jr — 10 4 - 5 GeV, and /_r = 10 5 GeV, respectively. Therefore, the R-axion leaves a displaced vertex 
from which a pair of the highly boosted muons is produced inside the beam-pipe in the first case, the inner detector 
in the second case, and the calorimeters in the last case. 



R-axion reconstruction 



We now consider the reconstruction of the R-axion. Since the two muons leave charged tracks inside the detector, 
it is possible to reconstruct the displaced vertex of the R-axion decay as well as the R-axion momentum. Notice that 
since the R-axion is highly boosted, the muons coming from the R-axion are collimated. The typical azimuthal angle 
between the muons and the R-axion is given by, 



J afi — • 
PT.a 



(45) 



where px, a is the transverse momentum of the R-axion. For m a = O(100)MeV and px = O(100)GeV, the angle 
<\>an is expected to be O(l) mrad. Therefore, the most important challenge to look for the R-axion is to separate the 
two muon tracks with very small open angle. Since a typical distance between two tracks at the pixel detectors is 
of O(100) /im, they can be separated by using the pixel detectors in principle. The actual separation efficiency and 
track resolution depend on the detail of the trigger system of detectors. In this study, we do not discuss the details of 
this problem, and instead, we simply use the track resolutions of a single muon track for each muon. At the ATLAS 
detector, for example, the resolution of the azimuthal angle of the muon is approximately given by, 

a ~ 0.075 8 - \ 8 _ (mrad), (46) 
Pt /GeVvsmfl 

where 9 is the polar angle of the R-axion. Thus, the azimuthal directions of the two muons can be well measured [41], 
and they are well separable. 

Related to the angle between the muons and the R-axion, the transverse impact parameters do of each muons are 
also expected to be 

d ~ — — CTuu <\> ai i — CT^. (47) 

m a 

The typical decay length at the rest frame of the R-axion given in Eq. (41) is of the order of 100 /im over a large 
parameter range. Therefore, the finite impact parameters of the two muons can be measured at the ATLAS detector 
where the resolution of the transverse impact parameter is given by [41], 

73 

a do ~ 11 ® (/mi). (48) 

PT/GeVysmo 
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Once the transverse impact parameter is measured, the transverse decay length is then determined by, 

L = d ° 

sin(0 ajU ) ' 



(49) 



This is expected to be several cm as shown in Fig.l. From the resolutions of the impact parameter and the azimuthal 
angle, the resolution of the transverse decay length is estimated by, 

^ = ^L ^. (50) 
L>t d (p a ^ 

Therefore, for d = 0(100) ^m and <p atl — 0(l)mrad, we expect the resolution of the transverse decay length to be 
marginal to better than 10%. 

Before closing this section, we mention the R-axion mass measurement. In terms of the observed momentum, the 
mass of the R-axion is given by, 

m 2 a -Ami +PnPP-°ln- ( 51 ) 

where p^) denote the size of the momentum and the transverse momentum of the muon (anti-muon). The absolute 
angle between two muons , 6^, is given by, 



A9 2 + sin 2 6 (52) 

where A9 and A<j) denote the differences of the polar and the azimuthal angles of the muon and anti-muon, while 6 
is the average of the muon polar angles. Since the R-axion is highly boosted, the expected value of 6^ is given by, 

6 ^~2sm6—. (53) 

PT,a 

At the ATLAS detector, by using the resolution of the azimuthal angle in Eq. (46) and the one of the polar angle [41], 

2 x 10~ 3 

<7 cot9 ^0.70x 10" 3 © — - , (54) 

PT/deV vsmv 

the resolution of the angle between muons is given by, 

^ ~ sm 2 9^-a cote © sin 2 0^^ © 

~ sin0^-cr cote © <?<$> © -> £), (55) 
2m a 2m a 

where we have used Eq. (52) and Ad = 2m a /p a , and A<p = 2m a /pT, a - Thus, the resolution of the absolute angle is 
expected to be of O(10) %. Therefore, the R-axion mass measurement is possible, although it requires careful study. 



V. BACKGROUND ESTIMATION 



In this section, we consider relevant background for the R-axion + jets events. As we have discussed, we select 
events with a jet of pr > 100 GeV balanced by a low-mass muon pair emerging from a displaced vertex. First of all, a 
potential background can come from mis-measuring the impact parameters of a prompt muon pair. The probability 
depends on the decay length of the axion and the tracking precision. Since we require the transverse impact parameters 
of both the muons to be around do in Eq. (47), the mis- measurement probability Pd is roughly estimated by, 

P ^ Erfc2 (7fc)' <56) 

where <Jd Q is the resolution of the transverse impact parameter given in Eq. (48). Thus, the probability Pd Q is quite 
small for the typical transverse impact parameter of the R-axion, do = O(100) fim, i.e. Pd a <C 10~ 5 . The limit of the 
suppression factor corresponds to do ~ 45^m, about 3cr from primary interaction vertex. 

The production cross section of a prompt muon pair with pt > 50 GeV is generically less than 10 pb. For example, 
the light unflavored mesons, such as p, w, <j), r\ can promply decay to two muons. The cross section of these meson with 
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TABLE II: Summary of rough estimation of muon pair background from different heavy meson decay modes. Here we have 
use mx = (5,2) GeV and cr = (460, 123) /im for B and D. As our signal cross section is about 10 fb, the background will be 
claimed insignificant if it's cross section is less than 0.1 fb. The final cross sections from the B° — + D ± + /i T + v — > + X 

and D° — » u> + K% — » nn + K% modes are further suppressed than given here (see discussion in the text). 

Pt > 100 GeV is similar to that of two jet event which is roughly cr mcS on = 1 Mb and the two muon decay branching 
ratio is -Brp™ eson ' = O(10~ 5 ) when taking into account the muon isolation. The muon pair can also be produced via 
the Drell-Yan process. The production cross section of such events is also about O(10) pb for px > 100 GeV. Hence, 
the cross section of the background event is less than about 0.1 fb which is much smaller than that of the signal (see 
Fig. 5 at the end of this section). 

More serious background comes from events which have a non-prompt muon pair. The main candidates for the 
background are as follows; 

• A muon pair in heavy flavor (B, D) meson decays. 

• A muon pair in light flavor (K) meson decays. 

• A muon pair in gluon splitting followed by meson decay. 

• A photon conversion in the detector material to a muon pair. 

Among these, the most serious background comes from the muon pair in a cascade decay of the B meson. The muon 
pairs from a cascade decay of the K meson and photon conversions will also be important if we look for the displaced 
vertex outside of the beam-pipe. We will estimate these background in turn in the following subsection. 

A. A muon pair from heavy meson decay 

The muon pairs from the heavy meson decays are serious background, since the decay length of these mesons are 
similar to that of the R-axion, i.e. ct = O(100) yum. Among them, some serious contributions are as follows; 

• A muon pair in non-resonant decays. 

• A muon pair in cascade decays. 

• A fake muon contribution in semi-leptonic meson decay. 

Non-resonant decays are inclusive meson decays of the type (B, D) — > + X . The Standard Model prediction 

of the non-resonant inclusive B meson decay into a muon pair is highly suppressed. The most serious mode is, 
B° — > K*°(i + n~ . For the cascade decays, (B, D) first decays to some resonances and followed by the decay of one of 
the resonances. When the 7r and if-meson decay in flight or punch through to the muon chamber, they can fake the 
muon. Thus, we should also consider semi-leptonic decay of heavy flavor meson with one of the light hadron (n, K) 
faking the muon. 

The effective cross section of each of these potential background events from meson X = (B,D) can be described 

by 

X X P. L, 

< Brffi x P Lt x P gco x Pp x nb, (57) 
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where we have used the fact that (b, c) 



100 GeV jet cross sections are about the same and they are <tb 



10 nb 



with pt > 100 GeV [42]. A factor P^ f represents the muon isolation condition that require no significant hadronic 

activity in the detector. It is generally estimated to be around 10 _1 . A factor Brffi is the effective branching ratio of 
decay channel. As we have to also take into account of the fake muon, this branching ratio should also include 
semi-leptonic decay with one of the light meson faking the muon. In the case of fake muon, the fake rate P^/^.k) 
should be included in the effective branching ratio 



A rough estimation of the fake muon probability P^/^.k) is given by 

rax Lt 



P n/X = -Pmis-id x Br x ^n+ 



v X [ 

Jo 



dL T exp 

ct x Pt.x 



Pt.x ct x 



xn x (L T ), (X = w ± ,K ± ), 



(58) 



(59) 



where P m ; s _id is the probability of the misidentification of the mesons to the muon, tix(Lt) is the rate of the punch- 
through mesons reaching to the transverse length Lt, and the r out is the outer radius of the calorimeter. In the case 
of the ATLAS detector, the calorimeter extends from an inner radius about 2 m to and outer radius 4 m, and the total 
thickness of the calorimeter is about 11 interaction lengths. By using these parameters, the punch-through rate at a 
given length Lt is given by, 



n x {L T ) 



e(L 7 



■>1 



(fa- 



ll 



exp 
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(r in -2m, r out ~ 4 m, Z\r ca i = r out - r in ). (60) 



As a result, we obtain rough estimations of the muon fake rates, 

^ ±A± < 2 x 10- 4 , 
P^/k± <4x HT 4 , 



(61) 



for pt > 50 GeV. In the above expression, we have assumed P m i s -id = 50% for tt^ and Pmis-id = 10% for [43]. 

The suppression factor P Lt represents the difference of the expected transverse decay length of the R-axion and 
the mesons. It can be estimated as Pl t — Lt/Lx where Lt and Lx are the transverse decay lengths of R-axion and 
meson X respectively. The factor P gco parameterizes a rejection factor based on the difference of the event geometry 
of the X meson decay and the R-axion decay. The factor P geo is estimated in the following way. A typical absolute 
angle between the muons from X meson decay is roughly given by, 



9 x 

MM 



~ 2sin( 



,- m x 
Pt,x 



(62) 



where pt,x is the transverse momentum of the X meson and mx is the mass of the meson. Thus, by comparing the 
typical angle between the two muon from R-axion decay in Eq. (53), the geometrical suppression factor is roughly 
given by, 



geo 



e x 

MM 



Px m a 
Pa m x 



( Pt,x m a 

\PT : a m X 



( Wa 

\m x 



(63) 



In the final equality, we have also used pt,x — PT,a, which is required for the transverse momentum balance between 
the muon pair and the remaining jets. In the above expression, the two powers of {0^^/9^) out of three come from 
the requirement that the absolute angle of two muons should be as narrow as 9^^, while the remaining power comes 
from the requirement that the transverse impact parameters of muons for a given transverse decay length arc O^ct^^) 
given in Eq. (41). 

A factor P^ is a rejection factor obtained from vetoing a 6-jet which does not fake the muon pair. For example, the 
6-tagging efficiency can be about 90 %, while the misidentification rate of a gluon jet as a 6-jet is about 0.1 [41]. Thus, 
we can suppress a background event by a factor of 0.1, while keeping the efficiency of the signal event at about 90%. 

The basic analysis given above does not include the momentum balance that requires the muon pair to acquire large 
portion of the jet momentum in order to pass the pt cut and that the sum of momentum of the muon pair points back 
to the primary vertex. We will use this basic analysis to eliminate most of the background and to identify dominant 
background which can not be eliminate by this basic cut. Since our signal cross section is about 10 fb, we will claim 
the background insignificant if its cross section is less than 0.1 fb. For those backgrounds which are lager than 0.1 fb, 
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FIG. 2: The momentum fraction of the visible energy of a 6-jet. The solid line denotes the visible momentum, i.e. pt — p v 
before smearing, while the dashed line includes the smearing effect. In the figures, we have assumed the inclusive neutrino 
mode of the 6-jet is 20 %, and the jet momentum resolution about 6 % for pr — 100 GeV. 



we will treat them more carefully in the next step. For m a < 1 GeV, where the displaced vertex is Lt = 1 — 5 cm, we 
summarize the estimation of various backgrounds in Table II. 11 

As been illustrated in Table. II, a muon pair coming from a cascade decay of the B meson, B° — > Z? ± + /i T + v — > 
K Q +ii + /.i^ +vv is the dominant background. Note that, although the subsequent D meson decay also leaves a displaced 
vertex, the typical distance of the vertex from the muon track from the _B-decay is too short to be distinguished. For 
these background events, we need more careful analysis in order to distinguish them form the signal. For the cascade 
decay of D meson, the background is about 1 fb. However, the muon-pair isolation cut in a c-jet is expected to be 
more efficient than that in a 6-jet, since the D meson momentum fraction to the c-jet is typically much smaller than 
the one of the B meson to the 6-jet. Thus, although we do not pursue further analysis, we expect that the background 
cross section from D° — > w + Kg is more suppressed than that shown in Table. (II). Further suppression can also be 
achieved by more detailed analysis similar to that given below for the B decay. 

In addition to the suppression discussed above, the background from B meson cascade decay can be suppressed 
further by requiring transverse momentum balance between the muon pair and the remaining jets. In Fig. 3, we show 
the momentum fraction of the muon pair in the cascade decay, where we have assumed spherical decays of the B and 
D mesons. To obtain the x = p^^/pb distribution, we have used the ps/Pb distribution given in [47, 48]. The figure 
shows that the spectrum of the muon pair with p^^/pb > 0.7 is highly suppressed. Also, the other jet which does not 
involve the muon pair may decay into neutrinos. In Fig. 2, we have plotted the visible energy fraction of the 6-jet, 
Pvis = Pb — Pv hi the figure, we have assumed the inclusive neutrino mode of the B meson, B —> D + v + anything, 
is 20%. From those distributions, we can obtain a rough estimation of the size of the missing transverse momentum. 
As a result, we found that the background is rejected by a factor better than 10 -2 , by requiring that the missing 
transverse momentum is smaller than 20 % of the observed pt of the muon pair. On the other hand, the signal 
efficiency by the missing transverse momentum cut, j>T < 0.2px.^ M , is close to 1, since the main cause of the missing 
transverse momentum in the signal event is the resolution of the jets, <7j e t ~ 0.6 GeV x ^p T /GeV [41]. Thus, for a 
conservative estimation, we can use Pa ~ 10~ 2 . 

Putting all these suppression factors together, the transverse decay length distribution of the background events is 
roughly estimated by, 

%^ <BrW xP xJ>x Up T ^^e^ 

dL T W ' P Jp T , mzn CT B PT.B 

where Pa ~ 10~ 2 and P^ ~ 10 _1 . In Fig. 3, we show the resultant transverse decay length distribution for m a = 

300 MeV. The figure shows that the final suppression factor in Eq. (64) is about 10 -3 for Lt — 1 cm. In the analysis, we 
have used the production cross section of the 6-quark pair production given in Ref. [42] with pb = 0.8pb- By comparing 



rriB Lt 
Pt,b ct b 



m a 
m B 



dai 
dp 



T.B 



(64) 



1 In the table, we have not included the exclusive decay mode B — > whose branching ratio in the Standard Model is too small 

(Bri^ = O(10 — 9 ) [45]) to give a serious background. In the table, we have also used Pgeo ^ 10 — B for m a ^ 1 GeV in the J/t/} mode. 
This comes from the fact that the exclusive decay J/ip — > (J-fJ. is a two body decay, and hence, the narrow opening angle between the 
two muons in Eq. (45) is difficult to achieve. 
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B -> \x\X + X 
m B =5 . 3GeV 
cr B =0 . 045cm 
m a =3 0MeV 

p T , mlm = 100GeV 



p T , mim =200GeV 



W. 

0.2 0.4 0.6 0.8 1 0.01 0.1 1 10 100 1000 

X L T /cm 

FIG. 3: Left) The momentum fraction of a muon pair to the B meson momentum (solid line) and to the 6-quark momentum 
(dashed line) in a cascade decay B° — > D ± y^u — > K° p,p,vv. In this analysis, we have assumed spherical decay of B- and D 
mesons. Right) Transverse decay length distribution of the B meson decaying into a pair of muons which fakes the muons 
from the R-axion decay. The solid (dashed) line denotes the distribution with a lower cut on the transverse momentum of the 
R-axion candidate, pT,min = 100 (200) GeV. The distribution scales by (m a /300 MeV) 3 for a different R-axion mass. 

with the transverse decay length distribution of the R-axion in Fig.l, it is found that the signal-to-background ratio 
is good especially for Lt > 0(1) cm and m a < 1 GeV. 12 A simple numerical simulation of the Pu+u- distribution in 
b cascade decay using CalcHEP[49] is also carried out and used to estimate the rejection factor. By imposing the cut 
that the invariant mass of the two muons is less than 1 GeV and the pr of both muons is larger than 40 GeV, the 
total cross section of pp — > bb — > fi + pT is about 0.15 fb using Brffi — 10~ 2 . After including the factor P^ ~ 0.1, this 
background can be suppressed to less than 0.015 fb. 

B. A muon pair from K meson decay 

The K mesons have long lifetimes: ctk — 15.3m for K\ and ctk — 2.68cm for Kg. Their inclusive decays into 
a pair of (fake) muons contribute to the background, if they decay at a short distance. The dominant decay modes 
which contribute the above events are, K\ — ► ^ + /i T + v (Br = 27%), K s — > 7r ± + p T + v (Br = 4.7 x 10~ 4 ), and 
Kg — > 7r + +7r~ (Br = 68%). Since those mode requires ir — > fi misidentification, the effective branching ratio is given 
by 

7T + /i + Z/) ~ 10" 9 , 
Tl + H + V) ~ 10~ 9 , 

tt + tt) ~ 3 x lO" 10 , (65) 

for pk,t ^ 100 GeV and r p i po ~ 5 cm. The cross section of the background is also suppressed by the requirement of 
the momentum balance and the muon pair isolation. For example, if we are require pn/Pict > 0.7, the background is 
suppressed by 0(1O -2 ) (see the K meson spectrum in a light jet given in Ref. [50]). 

Therefore, the background events are suppressed by 10~' 11 ~ 12 - ) , and hence, the background from the K meson 
decay is subdominant for px > 100 GeV where the cross section of jets is less than about 1 /ib. As a result, we find 
that the background cross section from the K meson decay is also subdominant. 




Br^(Kl^^)+n + v) = ( mK r ™ e )xP^xBr(Kl 

\PT,K ctk l J 

Br^\K% - Mtt) + At + „) = ( mK rpipe ) x x Br(K° s 

\Pt,k ct Ks J 

Br^\K%^^)+^)) = (2!L?SiBl)xI* /w xBr(K° s 

\Pt,k ct Ks I p/ 



In this analysis, we have not used a muon pair isolation cut on the background events. This requirement is not independent of the 
momentum balance requirement. 
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C. Heavy flavor meson pair from gluon splitting 



When the gluon splits into a pair of b quarks, we may have two B mesons in a single jet. In this case, both B and 
B can decay to muons, and this muon pair can contribute to the background if the two tracks intersect at the finite 
transverse length. The suppression factors to the background cross section of the accidental muon pair are as follows. 

* Br^ +X , the branching ratio that the both the B^ meson decays include a muon. 

* f26/giuom the angle between B mesons is narrower than 2itib/pb so that the two muons would intersect. 

* P^ T , the suppression from the missing Et cut. 

* -Pgeoj the muons intersect with an angle between muons of 0(2m a /ps). 

The dominant contribution to the inclusive muon mode of B ± comes from the mode B^ — > D° + /j, + v + X which 
the branching ratio Br ~ lO^ 1 . Thus, the background cross section is suppressed by Br 2 = 1CP 2 . 
The second suppression factor is estimated by, 

cr 2 6 1 /" Q ™ ax dQ 2 ,~J 2m 2 \ I 4m 2 " 

where cr g i uon ~ 1 /xb is the production cross section of the two gluon jet for px > 100 GcV. Here, the upper limit of 
the integration, Q 2 Max , is given by, 

4L* = 4ml + \pI ( 2— V , (67) 
4 V Pb J 

which comes from the requirement that the two b-quarks are collimated in an angle 2rriB/pB- By using the masses 
tub = 5.3 GcV and = 4.2 GeV, we find fWgiuon is about 10~ 3 . In this expression, we have assumed that the gluon 
multiplicity in a jet is 1, otherwise the momentum fraction the B meson pair to the gluon jet gets smaller. If we take 
into account the gluon multiplicity, the cross section is further suppressed (see Ref. [51] for detailed discussion of the 
gluon multiplicity). 

The third suppression factor P^ is significant in this background. Similar to the previous background from the 

cascade decay, B — > D + fi + v, the momentum fractions of the muons are suppressed in the inclusive muon modes, 
since the dominant contributions to the inclusive muon mode involve the D mesons. For example, if we require that 
the distribution of the momentum fraction of the muon pair carrying at least 70 % of the B meson momentum, is 
less than 10 -2 for each muons with ps > 50 GeV. In addition, since the accompanying jet is a gluon jet, the neutrino 
emission rate in the remaining jets is much suppressed compared with the 6-jet case. Thus, we can expect that 
the rejection factor from the missing transverse momentum cut, i.e. fx < 0.2pT iWt , is much better than 10~ 3 , i.e. 
Pa < 10- 3 . 

The suppression factor P gco is much more complicated compared with the other background events discussed above. 
However, since the angle between the two muons must be of 0(2m a /p a ), while the the muons from the B mesons are 
emitted in directions with angles of 0(ms/pB), the factor P geo at least involves, 

1 / m„ 



Here the factor 1/3 2 comes from the requirement that the both of muons are emitted in a face-to-face direction so 
that they would intersect. 

Put the above suppression factors together, the cross section for the background events are roughly estimated by, 

0-2B^2 Al +X < Br 2 B ±^ ll±+x X P gco X X P 2 b/gluon x (7 gluon , (69) 

and the total rejection factor is expected to be much better than 10 -10 for m a < 1 GeV. Therefore, we expect that 
the background from the combinatorial muon pair from two B meson decays is subdominant. 

Similar to the B meson background, a combinatorial muon pair from two D mesons contribute to the background 
if they intersect. By replacing 6-quark and B meson masses with the c-quark and D meson masses in Eq. (66), we 
obtain a suppression factor, 

P2D/ S luon - 3 X 10- 3 . (70) 
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7 xl(T 8 


Si 


14 28 2.33 


4.98 


2 xl(T 7 



TABLE III: The probability for the muon photo-production at the 1 mm thick material. In the probability, we have used 

Cyfj,(_i = 3. 



Pg 0O should also be modified as 

This lead to a suppression factor of 10~ 9 which is a little bit weaker than that of the B meson. Hoever, the muon 
isolation cut is more effective compared to the B meson decay event. Therefore, again, we expect that the combinatorial 
background from the D mesons is also subdominant. 



D. Muon photo-production 

The photo-production of a muon pair at the detector material is also a serious background if we look for the R-axion 
decaying outside of the beam-pipe. The cross section of such events are given by, 

^7 — >[ifi By — >[i[i X (Jj 7 (72) 

where ct 7 is the production cross section of the photon, P-y^^ is a probability of the muon photo-production for a 
photon for a given material. 

The probability P 7 ^ w is roughly estimated as followings. The muon photo-production occurs by exchanging a 
photon in a i-channel between the muons and the target nucleus, and the cross-section is dominated by the process 
with the minimum momentum transfer of order of 0{rnj l /k) for a given photon momentum k. Thus, we can take the 
the nucleus of the material as a point-like particle and neglect the recoil of the nucleus. In this limit, the differential 
cross section is approximated by the Bethe-Heitler cross section (see Ref. [53] for more discussion), and the total cross 
section of this process is almost insensitive to the incoming photon momentum; 

Cphoto-production = Cj^^ X 10 3 °Cm 2 X (^~^J ■ (73) 

where Z is the charge number of the nucleus. The small momentum dependence is factorized in a coefficient C 7AlM = 
O(l). From this cross section, the photo-production probability per a incoming photon at the material of thickness 
Al is given by 

Bj^^^ Cphoto— production X 71 X (74) 

where n denotes the number density of the nucleus of the material. In table III, we list the photo-production probability 
at the 1mm thick of Be, C and Si which are typical material of the beam pipe and the inner detectors. From the 
table, we find that the probability of the muon photo-production at 1mm thick materials is O(10~ 7 ). 

Now, let us consider the photon production cross section. The one of the most important source of the photon at the 
LHC comes from the direct photon production via the QCD Compton process, qg — > q*y, and the quark annihilation 
process qq — > 57. In Fig. 4, we showed the production cross section of the direct photon at the LHC with a low px 
cutoff as a solid line. The figure shows that the direct photon production cross section is 0(1) nb for px > 100 GeV. 
In the same figure, we have also plotted the the muon photo-production cross section at the 1 mm thick material for 
Be, C, and Si as dashed lines. 

In order to compare the background process with the signature, we also plotted the cross section for the event 
with the R-axion decaying within the transverse length L max for a given low pr cutoff (solid line). The result for 
imax = 5 cm corresponds to the R-axion decaying within the beam pipe. In this case, the main backgrounds from 
the muon photo-production occur at the wall of the beam pipes which consist of Be. By comparing the figures, it is 
found that that the signal-to-background ratio is good if we chose pT,min = 100 GeV. 
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FIG. 4: Left) The production cross section of direct photon at the LHC with a low pr cutoff (solid line) at leading order. The 
dashed lines show the muon photo-production cross section at the 1 mm thick material, i.e. a x P, for Be, C, and Si. In this 
figure, we have used C Mfl = 3 again. Right) The cross section for the event with the R-axion decaying within the transverse 
length L max for a given low pr cutoff. The solid (dashed) lines show the cross section for the event in which the R-axion decays 
within 5 cm (50 cm) for fn = 10 4 ' 4,5 ' 5 GeV from the bottom up. 

For fn > 10 4 5 GeV, however, it is difficult to produce enough R-axion for L max = 5 cm even for the integrated 
luminosity 100fb _1 . In this case, we need to extend the region of the displaced vertex search. For example, if we 
choose L max = 50 cm, we expect a sizable number of the R-axion production. However, in this case, the amount of the 
material at which the muon photo-production occur gets much larger. Therefore, we need to see whether the muon 
pair candidate for the signal produced at the detector material or the in between of the material. For such study, we 
need more detailed detector simulation. 

Notice that the photon in jets also contributes to the muon photo-production cross section. In this case, the cross 
section for the background process is given by, 

(^jet—tfifj, = -fjct x Pf^ufi x "'jet j (^5) 

In the region for pr = 100—500 GeV, the cross section of the inclusive jets are about 7 x 10 2 of the direct photon cross 
section [54]. The jet rejection factor comes from the requirement of the momentum balance of the muon pair and the 
remaining jets as well as the muon pair isolation cut. For example, if we are require Pir/Pjct > 0.7, the background 
is suppressed by O(10 -2 ) (see Ref. [50]). Furthermore, the photo-production from the ir° is always accompanied by a 
photon which is not converted into a muon pair. Therefore, we can expect the jet rejection factor O(10 _ ( 2 ~ 3 )), and 
hence, the net contribution of the jets to the muon photo-production is comparable to the photo-production from the 
direct photon discussed above. 

Finally, we mention here the typical event shape of the muons produced by the photo-production. The typical 
invariant mass of the muon pair produced by the photo-production process is of 0{2m IJ ). Thus, it seems challenging 
to tell the signature from the muons from the photo production by using the invariant mass distribution of the muon 
pairs. 

E. Summary of the background estimation 

In this section, we gave rough estimations of the backgrounds to the R-axion signature. As we have shown, the 
most serious background comes from the inclusive decay of the B meson, B° — > _D ± + fi T + v — > + vv, if we 

look for the R-axion decaying inside the beam-pipe. In Fig. 5, we show a schematic figure of the detectability of the 
R-axion at the LHC experiments for a given values of m a and fu, with Cx fixed at the value Cl = 2. In the figure, the 
shaded regions "05 < 1 fb" and "050 < 1 fb" correspond to the parameter spaces where the integrated cross sections 
of the R-axion for Lt = 1 — 5 cm and Lt = 1 — 50 cm are smaller than 1 fb, respectively. In the both regions, we have 
assumed px > 100 GeV. Notice that the total cross section is almost independent of the R-axion mass. 

In the small fn and the heavy m a region, the lifetime of the R-axion is rather short, and the muon pairs from 
the light meson decays can be serious backgrounds. Even worse, the transverse impact parameter gets too small to 
be distinguished from the primary vertex in such region. We have shown the parameter space where the transverse 
impact parameter is too small to be distinguished as the shaded region on the bottom-right corner. 
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FIG. 5: Schematic figure of the detectability of the R-axion at the LHC experiments for a given values of m a and fa with 
C L = 2 fixed. 



In the figure, we also plotted the region where the signal-to-background ratio is worse than 10 for the background 
from the short-lived meson decays, by assuming that we are looking for the R-axion vertex with Lt = 1 — 5 cm and 
Pt > 100 GeV. In this case the most serious background is the one from B° — > + fi T + v — > fxfi + vv. As we see 
from the figure, this background is more serious for the heavier R-axion, since the suppression factor P geo in Eq. (63) 
is less effective for the heavier R-axion. 

Notice that the backgrounds from the B meson decays can be negligible for the longer decay length regions, for 
example, for Lt > 5 cm. However, in such longer decay length region, the background events from the K meson 
decay as well as the muon photo-production process at the detector material are significant. Thus, the R-axion search 
outside of the beam-pipe requires careful study. 

Finally, we comment on the R-axion search for the lighter mass. When the R-axion mass is lighter than the two 
times of the muon mass, the main decay mode of the R-axion is the electron pair mode. In this case, the R-axion 
mainly decay outside the detector, and the R-axion is virtually stable inside the detectors of the LHC experiments. 
Thus, the R-axion events just look like jets + missing pr events. The detection of the R-axion in such a mass region 
is quite challenging, since the background cross section from the jets + Z° — > vv is much larger than the R-axion 
production cross section. 



VI. CONCLUSIONS 



In this paper, we have discussed the detectability of the R-axion. The properties of the R-axion are interrelated with 
the nature of the supersymmetry breaking sector, the messenger sector, and the solution of the fi and B[i problems. 
Thus, the R-axion can be a powerful probe of the physics behind the SSM. In the minimal R-breaking scenario, the 
R-symmetry breaking scale can be as low as Jr ~ 10 4 ~ 5 GeV, and the R-axion has a mass of 0(100) MeV. As we 
have shown, in such cases, a sizable number of R-axions are produced at the LHC experiments and can be detected 
by searching for the displaced vertex left by the R-axion decay. As a result, we found that the we can detect the 
R-axion at the LHC experiments for m a ~ 200 - 1000 MeV and f R ~ 10 4 GeV (Fig. 5). 

In this study, we have assumed that the /z-term has a vanishing R-charge. When we consider models with a yu-term 
which has a non-vanishing R-charge, the mixing angles of the R-axion and the CP-odd Higgs bosons are different 
from the ones given in Sec. III. Especially, when the R-charge of the /i term corresponds to 2, the mixing through the 
B[i term vanishes, and the R-axion does not couple to the Standard Model fermions directly. In that case, the main 
decay mode of the R-axion is into a pair of photons for m a < 1 GeV and into a pair of gluons for m a > 1 GeV. In 
those cases, we need different strategies to search the R-axion. 

Several comments are in order. 1) In our analysis, we have only considered the R-axion with m a = O(100) MeV. For 
heavier R-axion, however, this particle can also decay into a pair of the heavier SM fermions. In that case, the R-axion 
does not leave a displaced vertex, and hence, we need another strategy to find the R-axion. 2) The R-axion searches 
in the rare decays of such as T(IS') and J/tp particles are also interesting possibilities. So far, the constraints on fn 
from those rare decay modes are given fa > 10 3 GeV coming from Br(T — > 7 + a) < 10 - ' 5 ~ 6 ) for m a < 1 GeV [55] 
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(see also Ref. [56] for a review of the experimental constraints on the axion-like particle). 13 3) Since the R-axion 
couples to heavy fermions such as top-quark and the gauginos rather strongly, the production process accompanied 
by those heavy fermions might be interesting channels. 4) At the LEP experiments, the dominant production cross 
section comes from e+e~ — ► Z{"f) — ► Z(j) + a via the anomaly coupling given in Eq. (34) and e+e~ — > Z — > h(H) + a 
via the mixing between R-axion and CP-odd Higgs bosons. For Jr > 10 4 GeV, we found that the cross sections of 
those processes are below fb, and hence, we had no chance to have produced the R-axion at the LEP experiments. 
5) At the Tevatron experiments, we need to search for a low px R-axion for the integrated luminosity 10 fb -1 (see 
Fig. 1), which requires careful study. 14 
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